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Abstract 

Belavin's Z„-symmetric elliptic model with boundary reflection is considered on the basis of 
the boundary CTM bootstrap. We flnd non-diagonal A'-matrices for n > 2 that satisfy the reflec- 
tion equation (boundary Yang-Baxter equation), and also find non-diagonal Boltzmann weights for 
the ^^'i J -face model even for n ^ 2. We derive difference equations of the quantum Knizhnik- 
Zamolodchikov type for correlation functions of the boundary model. The boundary spontaneous 
polarization is obtained by solving the simplest difference equations. The resulting quantity is the 
square of the spontaneous polarization for the bulk Z„-symmetric model, up to a phase factor. 



1 Introduction 

Integrable models with a boundary have been studied in massive quantum theories [Q, |[ |[ |4[ ||, |^ and 
half infinite lattice models ^, |l^ |Tl|, ^ |l^ . The boundary interaction is specified by the boundary 
S'- matrix for massive quantum theories [Q, and by the reflection matrix K for lattice models 0. The 
integrability in the presence of reflecting boundary is ensured by the reflection equation (boundary 
Yang-Baxter equation) in addition to the Yang-Baxter equation for bulk (i.e., without boundary) 
theory @. 

It was shown in Q that the boundary vacuum of boundary integrable theories can be expressed in 
terms of the vacuum and the creation operators in the bulk theory. In the explicit bosonic formulae 
of the boundary vacuum of the boundary XXZ model were obtained by using the bosonization of the 
vertex operators associated with the bulk XXZ model [|l5j . 

The quantum Knizhnik-Zamolodchikov equations [p^ pT| are satisfied by both correlation functions 
and form factors for bulk field theories ||l^ and for bulk lattice models ^ with the affine quantum 
group symmetry. It is shown in Q that correlation functions and form factors in semi- infinite XXZ /XYZ 
spin chains with integrable boundary conditions satisfy the boundary analogue of the quantum Knizhnik- 
Zamolodchikov equation ^ . 
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In this paper we study Belavin's Z„-syinnietric vertex model with integrable boundary condition, 
the boundary Belavin model. The i?-matrix of Belavin's model is expressed in terms of elliptic functions 
of the spectral parameter z so that the i?-matrix has doubly quasi periodicity. Thus we expect that 
the iiT-matrix of the boundary Belavin model also possesses appropriate transformation properties with 
respect to z compatible to those of the i?-matrix. We shall show that under such assumption the K- 
matrix of the boundary Belavin model is inevitably non-diagonal for n > 2. Our solution is diagonal 
for n = 2 but different from the one used in . 

On the basis of boundary CTM bootstrap |l^, we find that the correlation functions for the 
boundary Belavin model satisfy a set of difference equations, the boundary analogue of the quantum 
Knizhnik-Zamolodchikov equation. Furthermore, by solving the simplest difference equations, we obtain 
the boundary spontaneous polarization which turns out to be the square of that for the bulk Z„- 
symmetric model p2| . 

The rest of this paper is organized as follows. In section 2 we review Belavin's Z„-symmetric model, 
thereby fixing our notations. In section 3 we give two non-diagonal solutions to the reflection equation, 
one is a constant iiT-matrix, and the other is an elliptic if -matrix. Furthermore, we consider the 
boundary analogue of the vertex-face correspondence to discuss the connection between our if-matrix 
and the boundary weights of the aII}_i model In section 4 we construct lattice realization of the 
boundary vacuum states and vertex operators from the boundary CTM bootstrap approach. In section 
5 we derive difference equations for A'^-point functions of the boundary Belavin model. We solve the 
simplest difference equations with iV — 1 for free boundary condition to obtain the explicit expression 
of the boundary spontaneous polarization. The result gives the higher rank generalization of that for 
the boundary eight vertex model In section 6 we summarize the results obtained in this paper, and 
give some concluding remarks. 



2 Belavin's vertex model and the reflection equation 
2.1 Elliptic theta functions 

For a complex number t in the upper half-plane, let A^. :— Z + Zt be the lattice generated by 1 and 
T, and Et- := C/A,- the complex torus which can be identified with an elliptic curve. For a,b G M, 
introduce the Jacobi theta function 



a 
b 



(z, t) := ^ exp {7rV^(m + a) [{m + a)T -(- 2(z -f 6)]} 



(2.1) 



mG2 



Hereafter a positive integer n > 2 is fixed and we will use the following compact symbols 



a2/n + 1/2 
ai/n + 1/2 



(z,r), e(^\z) = ^ 



1/2 - j/n 
1/2 



(z,nT), 



(2.2) 
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for CK = (ai, ^1^®!^ and for j G Z^; and 



n— 1 n— 1 

The superscript (n) and the subscript n wih be often suppressed when we have no fear of confusion. 
The elhptic theta functions are expressed in terms of the product series 



/'j-2n. J-2n^3 I -tin . -i1n\'i 



(t2;i2)3 



where 



(t2;t2) 
OO oo 

(a;(ji,--- ,gfe)oo := J| ••• J]^ (1 - agl" 



(2.3) 



mi— 771;^— 



2.2 Belavin's vertex model 



Let V 



and {vi}i^Zn be the standard orthonormal basis of V with the inner product (vj, Vk) 



Let Vz be a copy of V with a spectral parameter z. The Z„-Baxter model is a vertex model on a 
two-dimensional square lattice £ such that the state variables take on values of Z„-spin. Each oriented 
line of C carries a spectral parameter varying from line to line. We assign a Z„-valued local state on 
each edge. Let 

k 





7-1 




Zl 



I 

be a local Boltzmann weight for a single vertex with bond states fc, I £ Z„. Arrows denotes orienta- 
tions of lines. We now define the linear map on t/ji ® called the _R-matrix as follows: 



i?^^1.^^2(t;^. ,g,„;) = J2 {v^(E)Vk)R{zi- Z2)fl. 



Belavin [ pT| considered the Z„-symmetric model satisfying 

(i) RizYji = 0, unless i + k = j + I, mod n, 

(ii) R{zy^j^piXl = for every i , j, k, I and p e Z„. 

In terms of two linear map in V 



gvi = uj Vi, 



hvi = Vi-i, 



where uj — exp{2TT^/ —1/n) , the conditions {2A) can be rephrased as follows: 



R{z){g®g) 
R{z){h(S)h) 



{g(^g)R{z), 
{h®h)R{z). 



(2.4) 



(2.5) 



(2.6) 



3 



Thus the i?-matrix of Belavin's Z„-symnietric model is of the form 

R{z)^^R{z), Riz)= J2 ^Uz)Ic®Ioc^- (2.7) 

Here Gn — Z„ <^ Z„, and la = g"^h"^ for ct = (01,02)- The normaHzation factor k(z) will be given 
lator. The coefficient function Uaiz) is determined by imposing the i?-matrix satisfies the Yang-Baxter 
equation 

^12(^1 - 22)^13(^1 - Z3)R23{Z2 - Z3) = i?23(z2 - 23)^13(^1 - Z3)Rl2{zi - ^2), (2.8) 

where Rij{z) denotes the matrix on V®^ , which acts as R{z) on the i-th and j-th components and as 
identity on the other one. Belavin's solution to (|2.8|) is given as follows: 



Ua{z)^u^a\zM :^1 M2 + Wn) ^ 
n act\w/n) 

where 'w{^ mod A^) is a constant. It is obvious that the following initial condition holds: 

i?(0)=P, P{x®y)=y(^x. (2.10) 

In order to facilitate the derivation of the similar results for the if -matrix of the boundary Z„- 
symmetric model, we give brief sketches of proofs of several well known properties for Belavin's i?-matrix. 

Proposition 2.1 The Boltzmann weights or the elements of R-matrix are given as follows ^ 

h{z)e(^-'^^z + w) 

., , —7- — , , , , — 77 — T- — - if I + k ^ 1 + I, mod n, 

R{z)fi^^ 90-k){z)9('-3){w) ^ J , , ^2.11) 

otherwise. 

[Proof] Because of the Z„-symmetry, 

i^i^-'{z)^Kli{z)= J2 ^c^i^WodWc'ti = s'^+i E u^o.,,-,){^W^-'^"'. 



Set 7^"''(z) = Roa+bi^)- Then we have 

nt{z,v)= J2 <l,-a)(^>^V"'"^- (2-12) 



The transformation property of Tl°'^{z) and the initial condition TZ°^^{Q) ~ 5^^ imply that 



TZ'^^{z) = 0, at z = CT (c 7^ —6, mod n) and z = (a — h)T — w, mod A„,-. (2.13) 
Hence TZ°''^{z) has the form 

n''\z) ^ c''\w)e''''-''\z + w) W e^^Xz). 
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By substituting z — —br we have 



which concludes that ( p.ll ) holds. □ 



As a corollary of Proposition 2J wc have [p4[ 

PR{-w) = i?HP = 



(2.14) 



Now we assume that 0<t<q<u<l, where t := exp(7r-\/— Ir), q :— exp(7r\/— Iw), and u := 
exp(— 7r\/— Iz). Following Baxter p4| we call such domain of parameters the principal regime. Note 



that ( 2.11 ) is weights of the eight- vertex model when n = 2. 



2.3 Unitarity and crossing symmetry 

Belavin's i?- matrix satisfies the unitarity and crossing symmetry relations HJ, |2^, [2^ . 

Proposition 2.2 Belavin's R-matrix satisfies the following unitarity relation or the first inversion 
relation: 



R2l{z)Rl2(~z) = pi{z,w)I ®I, 



whe 



pi{z,w) 



a{z + w)a{—z + w) 



(2.15) 



(2.16) 



[Proof] Note that 



R2l{z)Rl2{-z) = ui^''(z,w)I-^ U^^\-Z,w)Ili®I^'' 



aeG„ 



= Y u^^\z)u^-\-z)I^^I(,®I^r^^ 



af3eG„ 



a£G„ 



where 



f^a\z,w)= J2 ^^'^'"'^u^£\z,w)u^:U-z,w), 

aeG„ 



and {a, a) — aia2 — Oi2ai. Proposition p.2| is thus reduced to 

fa\z,w) = pi{z,w)5aO- 



Concerning the proof of (2.18), see Theorem 3.3 and Lemma 3.2 in |2q|. □ 



(2.17) 



(2.18) 
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Next we describe the crossing symmetry for Belavin's Z„-symmctric model. For that purpose let us 
recall the i?-matrix on K (g) L, where K = Vz-^ (8) • • • ® Vz^. and L — Vz' ^ ■ ■ ■ ^ Vz' 



R 



V, V I V , 

-^l;fe+l ' 

TdK,L ID I JD 



'•l---fc;fe+l- 



YBE holds for R^'^ by virtue of YBE for R^'^ (2.8) 



rfK,L pK.M pL.M pL,M rfK.M tdK,L 
^12 ^13 ^23 ~ -"-23 ^13 ^12 ' 



(2.19) 



as a linear map on K ^ L (g) M 
For special = Vz^ g) ■ ■ ■ i 
associated with is given as follows |2 



such that zj = z + {k + 1 — j)w (1 ^ j ^ fc), the fusion operator tt 



^ ■— -^k-l:k -"•fc-2.fe-l:fe ■ ' ■ ^l,--- ,fc-l;fc 



From the first equation of ( ^.14 ) and the Yang-Baxter equation (2^) we have 



(2.20) 



(2.21) 



Let V* be the dual space of V and {v*}i^z„ be the dual basis of {vi}i^Zn- Then we have the 
isomorphism C : V*^^^^^ — > Anti(i4r"^^) 



^ 7(^7^ 



(2.22) 



where e^^ is the n-th order completely antisymmetric tensor. The spectral parameter z+nw/2 

associated with the dual space V* refers to the mean value of n—1 spectral parameters z+{n — l)w, • • • , 
z + w of V^. Then the i?-matrices onV^V* and V* i^V are defined as follows: 

JlV.,,V:^+„^/2 = (/(g,C')-li?^-l'^-2 + ('.-l)™®-®^-2+»(/® c), 

RK+^^/2y-2 ^ (c® /)-ii?^-i+<'-i)"'®-®^-i+™'^-2(C(g)/). 



(2.23) 



The un-normalized R on V V* and V* 'E)V are also defined in a similar manner. 



Proposition 2.3 The R-matrix on V®V* andV* ®V defined in ( 2. 2S ) meet the crossing symmetry 



p^.2.V;,+„„/2 (-n^^iy- ^ Vr h{-Zi+ Z2+pw) 

■^12 



-i^,2,y,i+„„ TT K-ZI + Z2 -pw) 

"^^ -'■ , h{w) 

p—i ^ ' 



(2.24) 



where ti denotes the transposition of the i-th space. 

^Notc that the spectral parameter of V* is shifted by nw/2 from the one in p^ . 
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[Proof] Let 

i,k 

Because of the initial condition ( 2.1C ) and the second equation of ( ^.14 ), the clement (— z) vanishes 
at — z = pw, where p = 2, • • • , n — 1. Thus we have an entire function from a*^'(— z) devided by 

h{-z - 2w) • • • /(-z - {n - l)w). 

The transformation property of z) are the same as Rj^j{z). It follows from the second equation 



of ( [2.14 ) that 6*1 (—z) — &t z — CT ior c ^ i — k and at z = (z — fc)r — w, which coincide the zeros 
of i?fcj (z) (|n|). Thus 6;|(-2) equals up to a scalar factor, which is determined by substituting 



z = (fc — i)T. The second equation of (2.24) can be shown in a similar way. □ 
From ( 2.16D and (2.24), we have the following second inversion relation ||2^ 



E^12(^)i?2\(- 



z — nw) = P2{z, w)I, 



where 



h{—z)h{z + nw) 
KHw) ■ 



(2.25) 



(2.26) 



Imposing the unitarity and crossing symmetry condition with respect to the normalized i?-matrix: 



(2.27) 



R21 



(<-'^-)*S 



R12 



(■^21 



(2.28) 



the normalization factor k(z) should obey the following functional equations: 



k(z)k(— z) 
k{z)k{—z — nw) 



pi(z, w), 
P2{z,w). 



Hereafter k{z) is often denoted by k(u) through the relation u = exp(— 7^^ 
regime using (2.3) the following expression solves ( |I29| ) H 



{u^;t^)o,{t^u-^;t^). 
{q';t^)oo{t^q~';t^)c 



-k{u), 



(2.29) 

Iz). In the principal 
(2.30) 



where 



k{u) 



((j2+2«y-2. ^2^ g2n)^(y2. ^2^ g2n)^ (^2^-2+2^^-2 . ^2 ^ ^2« )^ (^2^2 . ^2^ ^2n) 

From k(1) = 1 the initial condition for R also holds: 



R{Q) = P. 



(2.31) 
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3 Boundary Belavin model 

3.1 Reflection equation for the boundary Belavin model 

In this section we consider the following reflection equation or the boundary Yang-Baxter equation: 

K2{Z2)R21{Z1 + Z2)Ki{zi)Ri2{zi - Z2) = i?2l(2l - Z2)i^l (zi )i?12 (zi + Z2)K2{z2). (3.1) 

The reflection equation is valid when zi ~ Z2 because i?(0) — P. Furthermore, the following 

Lemma holds: 

Lemma 3.1 The reflection equation ( ^.l| ) is valid when (1) Zi ~ 0; (2) zi = —Z2 provided 

(1) Boundary initial condition: K{Q) — I\ 

(2) Boundary unitarity relation: K{z)K{—z) = /, 

respectively. 



[Proof] It is evident from the unitarity (2.27) and the initial condition (2.31) for i?-matrix. □ 



Here we notice that Belavin's i?-matrix have the following quasi-periodic properties 

R{z + l) = -~{g ® I)-^R{z){g ® I) = {I ® g)R{z){I ® g)-'^ , 

R{z + t) = -(/i(g)/)-ii?(z)(/i®/) X exp{-27rV^(z + f + f)} (3.3) 
= -{I®h)R{z){I(g)h)-^ X exp{-27r\/^(z + § + f)} . 

Thus we have the following Proposition: 
Proposition 3.2 Let 

K(z) = -f-ir(z), 



where \{z) is a scalor function. Suppose (SA) and the following quasi transformation property: 



K{z + l) = -gK{z)g, 

'K{z + t) = -h'K{z)h X exp{-2TTy^ [z + ^ + c)} , 



where c is a constant. Then K(z) solves (3.1) 



[Proof] Let F{zi, Z2) stand for the difference of the LHS and the RHS of ([3.1|). Then we have 

F{z^ + \,Z2) = -ig®I)F{zi,Z2){g®I), . 

, (3.5) 

F{zi+T,Z2) = -{h'Si I)F{zi,Z2){h'g) I) X exp{-2T:y/^B), 

where B — 'Szi + 'St/2 + 2w/n + c. The second equation of (^^) implies that the (ifc, j7)-th element of 
F(zi, Z2) satisfies 

F(zi + T, Z2)fi = -Fiz, + T, Z2y+fi X (-27^^/^i3). (3.6) 
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Thus wc find that F(pT, Z2)fi oc F(0, 22)-+^, = for ^ p ^ n — I from Lemma |3.l| . Similarly, wc have 
F{z2 + PT, Z2)fi = F{-Z2 + pr, Z2)fi = Ofor0^p^n-l: 

F{pT, Z2)fi = F{Z2 + pr, 22)}? = F(-Z2 + pr, Z2)}? =0, (0 ^ p ^ n - 1). (3.7) 

Assume that F{zi + r, 22)^7 is not identically zero. From Richey- Tracy's lemma (see section 3 in 
pjl or Lemma 2.4 in pq|) we conclude that F{zi + r, .22)!/ has 3rt zeros in Em- whose sum is equal to 



nc — 2w — 3n(n — l)r — (« + j)T. The contradiction to (3.7) imphes the claim of this Proposition. □ 



3.2 Solutions of the reflection equation 



Under the assumption of the quasi periodicity ( |3.4[ ) compatible to (3^) we find that the K{z) is not a 
diagonal matrix for n > 2. When n = 2 we can take K{z) diagonal because of g^^ = g and h^^ = h. 
The most general and non-diagonal solution for n = 2 is given in ^ . Other non-diagonal solutions 
for Z^i^-'-vertex model are given in 



In this paper we consider the following two solutions of (3.1), which can be also found in |31 



3.2.1 Constant if-matrix 
Proposition 3.3 Let 



ICoVj = Vn-i, 



(3., 



where Vn ~ vq- Then /Co solves (3.1). 

[Proof] It is easy to see glC^g — hICoh = JCq. Hence we have 

^2(2:2)^21(2:1 + Z2)Ki{zi)Ri2{zi - Z2) 
= / ® /Co ^ Uoc{zi + Z2){Ia^ ® Ia){ICo ®I)Y^ Uf^{zi - Z2)(/^ (g) Ip^) 

= /Co (g) /Co ^ t^"i"^Uc.(zi + Z2)Ia (E>IcJ2 ^^''^^Ufiizi - Z2)If3 ® /_/3 
OL (3 

= /Co®/Co^w''^^^M/3(2l -22)/^®/-/3^w"'"'Uc.(zi +Z2)/«(»/« 

= J2 '^^'"'M^l - Z2){I-(3 ® IfsKICo '^"'""^cizi + Z2){Ioc ® I-ocW ® /Co) 

(3 OL 

= R2l{zi - Z2)Ki{zi)Ri2{zi + Z2)K2{Z2), 

that implies this Proposition. □ 



3.2.2 Elliptic if-matrix 

Let 



TO 



n if n is odd, 
n/2 if n is even, 
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and let 

/c(z)= ^'"^"=4"c^(^,«)/2c.- (3.9) 

where 



for a = (ai, 0:2), and v(=/= mod A^) is a constant. Using the identity 

_j m— 1 

m ^ 



ai=0 



(5Q2_i if n is odd, 

<5Q2,i + Sa2,i-m if ?i is even, 



we have /C(0) = /Cq. 

Lemma 3.4 The following quasi transformation property holds: 
IC{z + l) = -.g-i/C(z).g; 

/C(z + t) = -/i-i/C(z)/i X exp{-27rV^(z+ § + ^)} 



(3.10) 



[Proof] This is based on the transformation properties of the ehiptic theta function. □ 
Lemma 3.5 Let K{z) = K.qK.{z). Then the boundary inversion relation holds: 

K{z)K{-z) = pi{z, v)L (3.11) 



[Proof] Direct calculation shows 
K{z)K{-z) 



where 



^0 Y 4a(^:^^)^a/Co Y ^^2^i-^^'")Jf3 

Y 4a(2>^^)^-a Y 4l3("^'^)'^/3 

Y 9a\z,v)Ia, 
a£G^ 



aeG„ 



(3.12) 



By comparing gl^\z,v) with f^'{z,w) defined in (2.17), we easily have g^'{z,v) — fa^^'{z,v) and 
hence ( 3.11 ) holds for even n. Repeating the similar argument in Proposition 2^ we can also obtain 
( pll) for odd n. a 



An) 



in). 



r{m). 



Theorem 3.6 Let K{z) = K,nJC{z). Then K{z) solves the reflection equation (3.1). 
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[Proof] From Lemma 3.4 we find that K{z) satisfies (p.4[) with c = v/m. Since K{0) — ICq^ — /, the 



K{z) also satisfies the first equation of (3.2). It follows from Lemma that K{z) satisfies the second 



□ 



one of (3^). Thus K{z) is a solution to the reflection equation (3T) from Proposition 3_^. 
Remark. Our ii'-matrix for n = 2 is different from the one used in so that the readers should be 
careful to compare our results with those of n = 2. 



3.3 Matrix elements of iC-matrix 

In this subsection we calculate the (j, A:)-th element of K{z): 

K{z)vu = 

Note that 



-{j-k)ai 



When n is even, thanks to the sum over a2, K{z)-'f. = if j + A: is odd. By comparing (|2.12| ) we obtain 



— , , TZm ^ ' ^ {z,v) ifj + fcisevcn, 

1 if j + fc is odd, 



(3.13) 



for even 



Tin 



I 



(z, v) j — k is even, 
■ (2, v) j — k is odd. 



(3.14) 



for odd n. 

We are now in a position to determine the normalization factor A(z). The boundary inversion relation 



( p. Ill) implies 

A(z)A(-z) = pi{z,v). 
Furthermore, the boundary crossing symmetry holds for n — 2 



which implies that 



j',k' 



(3.15) 



(3.16) 



(3.17) 



Since V* ^ A"-i(F) for n > 2, the LHS of ( |3.16| ) for higher n should be replaced by the (j, fc)-th 
element of the dual i^-matrix. We wish to discuss this point again in section 4. 
Here we assume the following functional relation holds for n ^ 2: 



(3.18) 
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Not ( 3.18| ) but ( |3.15 ) is important to calculate the spontaneous polarization in section 5, so that we 



proceed further under the assumption ( 3.18| ). By solving (3.15) and (3.17) we obtain 

1 ir^u^;e,q^-Uit\-^u^;t^q^^U ^i^' 



\{z) = 



where r = exp(— tt^— Iw), and 

(g"a;;t2,g2")^(tVa;;i',g'")oo 



(3.19) 



</.(x) 



(g2"j:; g2n -1^(^22.. ^2^ ^2«)^ (-^2^n2,. ^2^ ^2«)^ (^2^-2^™^. ^2^ g2n) 
(q^^+^X^; t\ g4n)^(i2g2n-2^2. ^2^ ^4„)^ 



(g2na;2; 54n)^(i2^2n^2. ^4„)^ ' 

3.4 Comments on boundary weights for the boundary A^^\ face model 

In this subsection we wish to discuss the boundary analogue of the vertex-face correspondence. Con- 



cerning the case n — 2, see [|13, 33 . Let us consider the bulk A^'^^-^-ia.ce model whose local state takes on 
values of P, the weight lattice of A^^-i I 
for a certain j G Z„, where 



An ordered pair (a, 6) G is called admissible if 6 = a + j, 

n — 1 



Vk. 



k=0 



Let 



W 



a b 
d c 



Z2 

a i b 



Zl- Z2 \ = 



--t--->~Zl 





b 












c 





be the local Boltzmann weight for a state configuration (a, b, c, d) round a face. Then W 

unless all the four pairs (a, 6), (a, d), {b, c) and (d, c) are admissible. Non-zero Boltzmann weights are 
given as follows: 





b 




(: 








c 





w 



where w{z, w) is a scalar function and 



W \ 

a + j a + 2j 



1 lab 

W 

w{z,w) I d c 



h{z + w) 



(3.20) 



W 



W 



a a + j 

a -I- j a + j + k 

a a + k 

a + j a + j + k 



h{w) ' 

h{ajkW — z) 



h{ajkw) 

h{z) h{ajkW + w 
h{w) h{ajkw) 



(3.21) 



ij + k). 
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Here 



a-jk = aj - ttk, aj = {a + p,Vj), 



and p — — 1 — j)j is the half sum of the positive roots. 

Jimbo, Miwa and Okado |Q introduced the intertwining vectors to show the equivalence between 
the Z„-symnietric model and the A^^\ model. Let 



(z + S ~ nwdj) if b = a + ij , 
otherwise, 



tf\z) := 



(3.22) 



where S is an arbitrary constant. Then we have the so-called vertex- face correspondence p2 

i?(zi-Z2)tS(zi)®<f(z2)=^M^( ^ Z1-Z2 I i^(zi)(»tg(z2). 



(3.23) 



Thanks to (3.23) the Boltzmann weights (3.21) solve the face-type Yang-Baxter equation 



b c 




( a 


b 




( f 


9 




Zl - 22 j W 






Zl - z}j W 


g d 






9 




{ » 


d 



fab 




' 9 c 




1 a g 




Z2 - 23 j W 


Zl - 23 j W 


\ 9 c 




\^ e d 




[ f e 



22 - 23 



21 - 22 



(3.24) 



Let us now consider the boundary A^^lj^-face model. By analogy with the bulk case, we find the 
following Proposition: 



Proposition 3.7 Assume that the existence of boundary weights V 's satisfying 









•) 







(3.25) 



Then V solves the face-type reflection equation 

a f 



22 w 



b e 



zi + Z2\V \ b 



Zl W 



a b 
d c 



Zl - 22 



(3.26) 



fa,e \ b g 



Zi^Z2\V\b^ 



lab 

zi\W , 

d e 



zi + Z2 \ V \ d ^ 



Z2 



In order to solve ( |3.25 ), let us recall the dual intertwining vectors jsj, |2^, |2[ 



a(r; 



t 



a{i) 



(2):= ($°(2))ydet$'^(2). 



(3.27) 
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Here is a matrix whose (i, j)-component is (2;), and $"(2) is a cofactor matrix of $'^(z). Note 

that ^^(z) is a column vector while t*J'{z) is a row vector. Thus by the rule of multiplication of matrices, 
f*''(z)ij^(z') represents a scalar function while t'^{z)t*J^{z') does a function-valued matrix. Since ij^(z) 
and ia''(z) enjoy the following orthogonal properties 



j=0 



the boundary analogue of the vertex-face correspondence (3.25) is equivalent to 

b 



V a 



z\ =c''(-^)i^(^)i?(^) = E^:o-)(-^)^Wi*c^'^W 



(3.28) 
(3.29) 



(3.30) 



Proposition 3.8 Let 



V a 



1—1 b 

V a 



X{z 



where A(z) is the same scalar function as for K{z), and V is defined by ( S.3(\ ). Then the boundary 
weights V 's satisfy the initial condition 



V a 



= <5,^ 



(3.31) 



and the inversion relation 



z\V \ a 



St 



(3.32) 



[Proof] The initial condition (3.31) follows from that for K{z) and ( 3.28| ). The inversion relation 
(I^H) follows from §J^) and ( ^ ). □ 

6 



The boundary weights V 



are non-diagonal in the sense that they do not vanish even for 



& 7^ c as a function of z. Hence ( 3.30| ) does not coincide with the diagonal solution of ( |3.26 ) involving 
the bulk Boltzmann weights for the A^j^2i-face model given in for n ^ 2. Such disagreement 



indicates that there may exist unknown solution to (3.1) corresponding to the solution given in [E3| and 



also unknown solution to ( 3.26 ) corresponding to our if- matrix, throughout the boundary vertex- face 
correspondence. 



3.5 Commuting transfer matrix 

The transfer matrix with L columns, 
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Vr V, 



1^1 



V2 Vx 



Z2 Z2 

is expressed in terms of i? and X-matrices as follows 



Z2 



Z2 



Tl{zi,Z2) = TroK+{zi)T{zi,Z2) 
T(zi,Z2) = TroT{-zi - Z2y^K^{zi)T{zi - Z2). 



Here 



T{zi - Z2) = R, 



T{-zi - Z2) 



= R 



01 

LO 



■R. 



■R 



e End (Fo ® Vi 
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"1 e End (Vb ® Vi 



Vl), 



are monodromy matrices satisfying 

-^12(21 - 2^2)71(21)7^(2:2) = 7^(z2)7'i(2;i)i?i2(zi - z[), 



(3.33) 



(3.34) 



and Tro signifies the trance on the auxiliary space associated with the spectral parameters zi and —zi. 
Note that the boundary monodromy matrix T(z, z') is a solution to the reflection equation: 

T2{z[,Z2)R2l{zi + Z2)Ti{zi, Z2)Rl2{zi - Z2) = R2l{z'i - Zl )Ti (zi , Z2 )-Rl2 (^1 + z[)T2{z[ , Z2) ■ (3.35) 



Proposition 3.9 // one takes 

K^{z) = K{z,v), K+{z) K{^z- ^w,v') e End{Vo), (3.36) 
where v and v' are arbitrary parameters, the transfer matrices ( 3.3t\ ) commute each other 

[TLiz,,Z2),TLiz[,Z2)]^0. (3.37) 



[Proof] From the crossing symmetry ( 2.28|) and the unitarity ( 2.27 ) we have 



Tl{zi, Z2)Tl{z[, Z2) 

= TriXi(-zi - fu.)ri(zi,Z2)Tr2if2(-4 - 'jw)T2{z[, Z2) 

= TtiTt2K2{-z[ - ^w)KlH-Zi - ^W)T^''{ZI,Z2)T2{Z[,Z2) 

= TtiTt2K2{-z[ - ^w)Kl^i-zi - ^w)Rl\{-zi - Z2- nw)R{\{zi + Z2)r^'{zi, Z2)T2{z[, Z2) 

= TriTr2if2(-4 " t^^)(i?2i (-^i - Z2-nw)K^{~z^ - ^^w))'^ {T^{z^, Z2)Ri2{zi + Z2)Y'r2{z[, Z2) 

= TriTr2i^2(-4 - tw)i?2i(-zi - Z2 - nw)Ki{-Zi - f w)i?i2(z2 - Zi) 

X i?2l(zi - Z2)Ti{zi, Z2)i?12(zi + Z2)T2{z[,Z2), 
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where we use TtAB = Tr Furthermore, from (p5| ) we have 



= TriTr2i?2i(22 - zi)Ki{-zi - jw)Ri2{-zi - Z2 - nw)K2{-z[ - |w) 

X 7^(2:^,2:2)^21(2:1 + 22)Ti(zi, Z2)i?12(2l - 22) 

= TriTr2i^i(-2i - '■^w){K2{-z[ - %w)Ri2{-zi - Z2 - nw)f- {R2i{zi + Z2)r2{z[, Z2)mzi, Z2) 

= TriTr2ifi(-2i - ^w){Ri2{-zi ^ Z2 - nw)K2{-z[ - %w)yHT2{z[, Z2)R2i{zi + Z2)Y^ri{zi, Z2) 

= TriTr2ifi(-2i - riw)Kl^{-z[ - ^w)T^^ {z[, Z2)r^{zi, Z2) 

= Tl{z[,Z2)Tl{zi,Z2), 



that imphes the commutativity (3.37). □ 



4 Boundary CTM bootstrap 

In this section we construct lattice reahzation of vertex operators and the boundary vacuum states for 
the boundary Belavin model. 

4.1 Partition function 

Let us consider the inhomogcneous lattice Clm with 2M horizontal lines carrying alternating spectral 
parameters 21 and —21 and L(= mod n) vertical lines carrying the spectral parameter 22 as below: 



The lattice Clm and the i-th ground state. 
The arrows stand for the orientation of the 
spectral parameters. The dots •'s stand for 
the boundary interaction K{z). 
For the sake of simplicity, we here denote 
the state i ± 1 and i ± 2 by i± and i±±, 
respectively. 

A zigzag line on which the state variables 
take i + 1 is presented for transparency. 

In this paper we restrict ourselves to the principal regime 0<t<q<r< u± < 1, where 
u± — exp (— 7r%/— 1(21 ± 22)). In this regime of parameters, the bulk Boltzmann weights of the type 
RizYj^j^l dominates the others; and the boundary Boltzmann weight Kl{z) is the largest among Kf{z) 
for fixed i. Thus in the low temperature limit t,q ^ Q, only the configuration such that the spin variables 
take the same value along the zigzag line (see the above figure) and increase by one in the direction from 
West to East, is possible. We call it a configuration of the ground state labeled by the boundary state 
i G Z„. Actually, the boundary weight Kq{z) (and K™{z) if n is even) are the largest among Kl{z). 
We therefore have only one real ground state for odd n and two for even n. Nevertheless, we regard all 
n kinds of configurations as the ground states. 
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In what follows, we fix one of them (say, i) and define all the correlation functions in terms of the 
low-temperature series expansion (i.e., the formal power series of t and q). Then the lowest order of them 
comes from the i-th ground state configuration. Furthermore, any finite order contribution is derived 
from the configurations which differ from that of the i-th ground state by altering a finite number of 
spins. It is equivalent to taking the GNS representation obtained from the i-th ground state (i-th GNS 
representation) as the Hilbert space. It is expected that the correlation function defined in such a way 
is an analytic function which has a finite convergence radius if there exists the phase transition at a 
finite temperature. 

Following we conjecture the partition function zj^*|^(zi, 22) of this model behaves in the thermo- 
dynamic limit L, M 00 as 



- LM(log^W(zi - Z2) +log^W(zi +Z2)) 
+ M (log (zi) + log (-Zl - f w)) . 



(4.1) 



Here ^^^\z) is the partition function per cite for the bulk theory, and v'^^\z) is the that per boundary 
cite, which are normalized as follows: 



«(z) = 1, v^°\z) = 1, z/(™)(z) = l,ff n is even. 



(4.2) 



Next we consider the boundary CTM lattice as below: 





The inhomogeneous CTM lattice split into 



four sections. 



r 
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and also denote the upper and lower lines of K{z) by i{B\ and \B)i, respectively. Let 

{•••«) Wp(3) «) -^^(2) «) Wp(i) bC?) e '^n,p{j) =i + l- j (mod n) for j > 1}, 



(4.3) 



:= {•••«) Wp(3) «) -^^(2) (g) Wp(i)|p(j) e 'Zn,p{j) = i (mod n) for j > 1}, 

and and 7i*(*) be their dual spaces. Then in the infinite lattice limit we conclude that \B)i e H, 

AB\ e n*^'\ and 



4'^(zi,Z2): 7YW~.hW, 
a(;V(zi,Z2): HW— .7i*W. 



The partition function is given as follows: 



Z^'^ (Zl , ^2) - . {B\A%^yy {Zl , Z2)^il (^1 ,^2)15). 



(4.4) 



(4.5) 



4.2 Vertex operators 

Let us introduce the type I vertex operators 



Zl- 



zi- 



Zl- 



Z\- 









i 








i 








j* 




^2 ^ 















= 0^._,_^)(zi-Z2):W«^H^'-i\ 



6;('~''')(z2-^i) 



Zl Z2 Z2 



where the sub/superscripts (i ± specify the spaces intertwined by the vertex operators. We often 
suppress these sub/superscripts when we have no fear of confusion. 

It follows from the Yang-Baxter equation that these vertex operators satisfy the following commu- 
tation relations |j: 

<j^Hz2W'{z^) = ^(i?^^l^^'=)j;j?0^'i(2l)(/.^"^(z2), 



<P*^^{z2W^{z,) = 5^(i?^-^-)j;j?0^i(zi)</.*^-^(z2), (4.6) 

r^-{z2W=^{zi) = ^(i?^^^-)^:;^|<^*^;(zi)<^*^n^2). 
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Furthermore, the unitarity relations for i?-niatrices imply the inversion relation of the vertex operators: 



(4.7) 



From the crossing symmetry we have 

r^iz) = (bji-z - ^w), ^ji-z) = 'P'iz - ^w). (4.8) 
Using these vertex operators, the transfer matrix for the semi-infinite lattice is defined as follow: 

Tb{zi,Z2) = (f>j{zi + Z2)Ki{zi)(j)^ [zi - Z2) 



j.k 

= ^ ^■'■(-^i - - Z2)Ki{z,)(b\z, - Z2). 

j.k 

If the i-th vauuam states jvac); and i(vac| satisfy the following reflection properties: 

^X^(z)0'=(z)|vac), = z.«(z)(/.J"(-z)|vac)„ 



(4.9) 



(vac|^0fc(z)4(z) 



j/(')(z)i(vac|(/)j(-z). 



(4.10) 



these vacuums are the eigenstates of Tb{z,0) associated with the eigenvalues ;y(*)(z), respectively: 

TB(z,0)|vac)^ = z/(^)(z)|vac)^, ,(vac|TB(z, 0) = i.(*)(z),(vac|. 

For n = 2, it is suffices to consider only two types vertex operators 4'-' [z) and 4>j{z) because of 
(jf^{z) — (j)i-j{—z — w) a,nd (l)*{z) ^ (j?-^^ {—z — w) Furthermore, from rB(zi, Z2) = 7s (—^i — Z2) 
for n = 2, we have 

<l>'^'i-zi - w - Z2)Ki{zi)cl)''{zi - Z2) 
= '^(f)^-''' {zi - Z2)K^! {-zi-w)(lx'' {-zi -w - Z2) ^^^^^ 
= J2 ^(-2^1 - w){2ik'^'''i-zi -yj~ Z2)cj>\zi - Z2)Kl,{-Zl - w), 



which implies the boundary crossing symmetry ( ^3.16| ). 

The crucial point in (4.11) consists in the self-duality (t)*j{z) = 0^~-'(z) for n = 2. Thus the boundary 
crossing symmetry ( |3.16 ) does not have a simple generalization for n > 2. We should rather regard the 
RHS of ( 3.1(]| ) for general n as the definition of the dual -matrix. In order to see that, let us repeat 
the reduction ( 4.11 ) for general n. Using eqs. (4^), ( 4.1C ), (4.6) and (|4.7|) we have 



i/«(0) 



''''\-z-^w)Kf,{z)cty>'{z)\Y^c)., 



= ^.(vac|. 

j' ,k' 

= ,(vac|0^(z)(i?^-^-— /2)/'=V|;'(z)</)*'=(-z - f u;)|vac). 



i,k 

j'.k' 



J2 .(vac|0*(-z - ^w){R''-''-^—f%''K^:{z)rH-z - IHIvac). 
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Thus, if we define the dual X-matrix by 



j',k' 



(4.12) 



then the following dual reflection properties hold: 

k 

,(vac|^0^(z)i^*(z)^^ = ^.W(-z-|u;),(vac|0*(-z). 



(4.13) 



The associativity condition of the algebra (4.6) and (4.13) implies the reflection equations involving 
iC*-matrices: 



K-({z,)B^;^'''-'K2{z2), 



V* V* V V' v v 



(4.14) 



4.3 Derivation of the reflection properties 

In this subsection we derive the reflection properties ( 4.10|4.1^ ). For that purpose we introduce following 
further two types of vertex operators: 



Zi- 
-Zl' 
Zi- 



-Zl' 

= ^^,_^^^)(zi,Z3) ^1- 

-zi- 

Z\- 



= iff ^^*'(zi,Z3) : — > 7Y("-i). 



zz 



where the sub/superscripts (i ± l,i) specify the spaces intertwined by the vertex operators. Hereafter 
we also suppress these sub/superscripts. 
From the reflection equation (13.11) 
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Zi- 



-Zl' 



z\- 



-Z\' 



z-i 

j 



\B) 



z\- 



-zi- 



Z\- 



^1 




-Z\' 



Zl- 



z/W(z3) X -^1- 



Zl- 



Z3 Z3 -Z3 

we have the following relation: 

k 

By similar argument we have 

k 

Furthermore, we have the relations 

4y)(zi,Z2)^^'(^l,^3)|i?). = cj,^{z3~Z2)Af^{zi,Z2)\B),, 

,{B\(pj{zl,Z3)A''^:^^\zl,Z2) = ,(B|A^\y(zi,Z2)0j(z2 -23)- 



(4.15) 



(4.16) 



(4.17) 
(4.18) 



These are based on the unitarity and Yang-Baxter relation of i?-matrix in the thermodynamic limit. 



The unitarity ( |2.28D allows us to obtain 



Z2 Z2 Z2 Z2 



\Bh 



Z2 Z2 Z2 Z2 



\B)i 



















-zi 










Z\ 






-z\ 






Z\ 





Z\ 



-zi 



Z\ 



zz 



zz 



Using the Yang-Baxter equation (^^) we get 
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Zl Z2 Z2 Z2 



-zi 



\B) 



Z2 Z2 Z2 Z2 



\B) 



Z\ 



-z\ 



Z\ 



z-i 



















-zi 










Zl 


i 
























Zl 



Z3 



By successive use of the YBE and the unitarity we can bring the hne associated with the spectral 
parameter Z3 to the directions pointed by dotted hues in the above figure as far as we hke. Thus wc find 



Z2 Z2 



Z2 Z2 



\B)i 



Z2 



Z2 Z2 



Z2 



\B)i 













3 

> 
> 










3 


















-Zl 




1 






-Zl 


> 








Zl 










Zl 








-Zl 




1 


-Zl 


> 




Zl 








Zl 



Z3 



Z3 



Those manipulation implies (4.17) because the contribution of Boltzmann weights along the tail graph- 



ically represented in the figure by the dotted line is unity in the thermodynamic limit. The relation 



(4.18) can be similarly obtained. 

Applying Ag^^ (^i, 22) (resp. ^^^^ (zi, Z2)) from the left (resp. right) to both sides of ( 4.15| ) (resp. 
(4.16)) and using (4.17) (resp. ( 4.1^ )) we obtain 



Y,Kiz,)i(l,'^{z3-Z2)A'^Uzi^Z2)\B), = iy^^\z,)<t>^~Z3-Z2)A^U^i,Z2)\B),, (4.19) 

k 

J2^iB\A^Nw(^l^Z2)Mz2 + Z3)Kizs)^ = u'^'^Z^UBlA^^^^izi, Z2)<Pjiz2 - z^). (4.20) 

fc 

Taking account of ( |4.19| ) and ( |4.20| ) with we find the following identification 

|vac), = Af^{zi,Z2 = 0)|B),;, ,(vac| = ,(B|4'^(zi, Z2 = 0). (4.21) 
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From the identification (4.21) and the definition of the dual if-niatrix (4.12) we obtain 

Y,K*{zs)ir'^{zs - z,)Af^{z^,Z2)\B), = ^.«(-Z3 - 'jwWi-zs ~ Z2)4w(^i. ^2)|i?).(4.22) 

k 

J2^{B\A'^Nwi^l,Z2mZ2 + Z3}K*{z3)^ - i/« (- Z3 " f U'). (S | ^j;)^ (^1 , ^2)</>* (^2 " ^^3 ) ■ (4.23) 



5 Correlation functions and difference equations 

The relations appeared in the previous section are not rigorous because all the objects are defined on 
the infinite lattice. Nevertheless we assume that eqs. (4.1- 4.2 are exactly correct on the basis of the 
CTM bootstrap method, which is supported by some numerical calculations [l4| and consistency with 
the vertex operator method p5[. 



5.1 Local state probabilities 

Let us consider the correlation function on the dislocated CTM lattice 



cy (z, z'\z[, ■ 



,z'^,ZN,--- ,zi)Ji'-J«.^".-.J"i 



-N 



Jn 

JN 



ZN 



Thanks to (4.17) and (4.18) we have 



z[ ^{B\ 



Ji 
ji 



Zl 



G%\z,z'\z[,--- ,z'^,ZN,--- ,zi)A^-'^'«'i-' 
,(B|4*]y(z,z')0,;(^' - z[) ■ --^f^iz' - z'^W«{zM - z') ■ • •0J"i(zi - z')Af^{z,z')\B) 



(5.1) 



23 



Thus the correlation function G^^ {z, z'|zj, • • • , z^, zjv, • • • , " ' '^^ normalized by the parti- 



tion function (4.5) is called the A^-point local state probability of the boundary Belavin model if we set 
zi ~ z[ = z' = 0, ji = j[ (1 ;^ Z ^ N). Owing to the unitarity we have 

Z«(zi,z2)= ,^jv,^iv,--- (5.2) 

Thus we obtain the expression of the n-point local state probability: 

5.2 Boundary analogue of the quantum Knizhnik— Zamolodchikov equation 

Only for n = 2, the A^-point function (^^) is reduced to the following 2A^-point function of the form 

F«(z,z'|yi,-- - ,yw,z^,--- ,zi)^^-.^«.^".-.J-i ^^^^ 
= ^{B\Af^iz, z')<j>'^^ (yi - z') • • • 0''^" (yAT - z')0^" (z^ - z') • • • ^-''^ (zi - zOA^!]^ (z, z')|i3)„ 

by putting yi = z[ — w and fc; = 1 — ][ for 1 ;^ ^ ^ ||^ . 

It is nothing to do with any local state probabilities for ti > 2, however, we can consider the 



correlation function of (5.4)-type: 

F^''(z,z'|zi,--- ,ZAr) = ^ Vj,®---®Vj^F'}i\z,z'\zi,---,ZNy^'-'^'' 



ji,---jN (5.5) 
'^'7z,z'|zi,... ,z^)^--'-'« = ,(B\Af^{z,z'W^iz,-z')---<jy^-izN-z')Af^{z,z')\B),. 



Here we assume that iV = mod n for simplicity. 
From the same discussion as in |l6[ ^ , we obtain 



Proposition 5.1 The correlation function ( |5.q j satisfies the following relations: 

1. R-matrix symmetry: 

2. Reflection property I : 

KNizN)Fl^\z,z'\zi,--- ,zn_i,zn) = iy'^^\zN)Fl^''{z,z'\zi,--- ,zn-i,-zn), (5.7) 

3. Reflection property II : 

Ki{zi)F^j^''{z,z'\zi,Z2r-- ,zn) = iy'^^\zi)F^j^'' {z, z'\ - zi - nw, Z2, ■ ■ ■ i zn), (5.8) 

where 

k{z)vk = Y,VjK*{-z-^w)]. 
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[Proof] The first equation (p^) follows from the commutation relation (4.6), while the second one (5.7) 



follows from (4.19). Finally, from the crossing relation ( pTq ) and (4.23) 
Ki{zi)Fj^\z,z'\zi,Z2,--- ,zn) 



E ^ • • • ® z')'/'*; -z,-^w)--- ^^^(z, z')\B),Kl{-z, - f 



^^^^H^i) E ® • • • ® z;,„.(B|4+^)(z, z')0;,(^' + ^1 + f z«) • • • 4w(^> ^')|i?>., 



we obtain the last equation ( |5.8| ). □ 

Owing to the equations (5^5.8) we obtain 

Theorem 5.2 T/ie correlation function ( [5. 4 ) satisfies the following difference equation: 



T]FI^\z,z'\zi,- ■ ■ ,zn) 



''-'■■■R.;' Kj{-zj) 



X i? 



3 + 1'^-., 



X i^,(z,)i?[^'''^" • • • R^;^:^:^'^^' F^^\z, z'\z,, • • • , z^v), 



where 



Tjf{z, z'\zi, ■ ■■ ,Zj,--- ,zn)^ f{z, z'|zi, • • • , Zj - nw, • • • , zn)- 



(5.9) 



Using the crossing symmetries we have another expression of the correlation function on the dislo- 
cated CTM lattice for general n ^ 2: 



, (B\J^^ (z, z')^*y^ (z* - z') • • • <jy*^'- {zh - z'W- (ZAT - z') • • • <P'' (^1 - ^O^sw l^> 



(5.10) 



where z* = z[ — for 1 ^l ^ N. We thus introduce the y*®" ® y®"-valued correlation function 

(z, z'jzj", • • • , Z^, ZAT, • • • , Zl) 

E 1"^^ '^t'jiv '^•••«'^'JlGj;^(z,z1zl^••• ,z;r,ZAr,--- ,zi)^i'---'^«'^" (5.11) 



Let us describe the i?-niatrix symmetry corresponding to (|5.6| ) . 
Proposition 5.3 Let 



= E 



y{N) 



,1N 



G^^'\z,z'\x^^,^,--- ,a;.(2iv))'"<^''-''"<^"' 
= .(i?|4^(z, z')<^>'^W . . . ^-(^Af^iz, z')\B), 

Here a be the permutation o/ (1, • • • , 2N), and 

z*^z'i-^w, (l^l^N); 
Z2N+1-U {N + l^l-^2N)- 



XI 



J2N+1-U {N + l-^l^2N)- 



(5.12) 



25 



and 



Then the following R-matrix symmetry holds: 

'* (^i ^'1 ■ ' ' J ^crU+l) J ^o-(j) !■■■)— ^cr{j),a{j + l) " ' ^cU) ' ^crU+l) J ' ' ' )i (5.13) 

where 

\v^, {N + \^l^ 2N); 
and Uj is the permutation of (I, ■ ■ ■ , 2A^) obtained from a by transposing a{j) and a{j + 1). 

The reflection properties can be similarly shown as before: 
Proposition 5.4 The following relations holds: 

K2Niz,)G^;;'\z,z'\--- ,z,) = iy^^Hz,)G^^'\z,z'\--- ,-z,), (5.14) 

k2N{zi)G^r^{z,z'\z,,---) = j.«(zi)riG^'^(z,zV^i,---), (5.15) 

kl{zl)G^^'\z,z'\zl,---) = ,,^^\~zl- ^w)T,G^^'\z,z'\- zl---), (5.16) 

K*,{zl)G^;'\z,z'\--- ,z[) = ,yi»){-,l-^^)Gp\z,,'\... (5.17) 



Here, 



^z-\w)). 



k*{z)vl^J2''*^K{- 

3 

and n, p, t £ &2N such that 

7r(2iV) = 2N, p{l) = 2N, <j(l) = 1, t{2N) = 1. 



[Proof] The relation (5.13) is evident from the commutation relations (4.6). The last two (5.16) and 



( p. 171) follow from (^.201) , (^) and ( [4.221 ). □ 
From Propositions ^.3| and 5.4, wc have 

Theorem 5.5 Let Vj = V-^,, 

, X2n) 



TiG''^\z,z'\xi,- 



for 1^1^ N, and 



TiG^^{z, z'\xi, ■ 



,X2n) 



forN+l^l^2N. 



Then the following difference equations holds 



Kt{xi)R^ 



■Rfr'''kt{-xi) 

^l-ll ^l+ll 



■R 



2Nl 



I ,,2N 



V',V 
'12N 



■ ■ '-^Tv (a^i, ■ • • ,X2n), 



ii-i 
11 



■ R^i^""' kii-xi) 



yl yl + l yl y2N 
' ' ' ^l + ll '''^2Nl 



k{xi)R^2N 



I ,r2N 



■R 



ii+i 



Gy{xi, ■ ■ ■ ,X2n), 



(5.18) 



(5.19) 



Theorem ^.5| gives an elliptic generalization of the corresponding difference equations for the bound- 
ary C/q(sZ„)-symmetric model [p2|. 
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5.3 Boundary spontaneous polarization 



Applying the similar argument as in ( |5.9D to the simplest case TV = 1 we obtain the following difference 
equations: 



TiG^(\z,z'\zl,Z2) = KI{-z*^)R;-^' -''^Kl{zl)R^'^^'''''G'^>[z,z'\zl,Z2) 
T^{z,z'\zl,Z2) - rII' 



''^ K2{-Z2)R,;^ ' ■ K2iz2)G'i^ (z, z'\zlz2) 



(5.20) 



where z* = zi — ^w. It is difficult to get each element of G^i \z, z'\zi, Z2), however, it is possible to 
obtain the expression of the following sums: 

n-l 

pW(z,z'|zi,Z2) = ^^™^Gf^(z,z'|zi - ^w,Z2)". (5.21) 
Note that the boundary spontaneous polarization as the vacuum expectation value of the operator g at 



boundary is expressed in terms of (5.21) as follows: 



is) 



_ PPiz,z'^0\z,,Z2) 



pI,'\z,z' ^0\zi,Z2) 



(5.22) 



Z1=Z2=Z' 



Now we restrict ourselves to the free boundary condition r — > 1 for simplicity. Since lim /C(0) 7^ /Co, the 
initial condition docs not hold if we take K{z) = JColC{z). Thus we should regard the if-matrix in this 
limit as K{z) = /C(0)/C(z). Under this identification the iiT- matrix behaves as 

K{z) k{z)In, 

where k{z) is a scalar function of z. 

Here we cite the following sum formula from ||2^ p| 



9'•^^z + w) h{{z — m)/n + w) 

— J~l 



Then we see the dual X-matrix in the free boundary limit r — s- 1 behaves as 



(5.23) 



K*{z - f w) ^ fc(-z)/o(u2(j")/„, 



where 



fM := ^c.™^"i?(z) 



1 (u;-"'g'u-"/";t^)oo(^'c^"^g-^M^/";t')o, 

The difference equations ( 5.2C ) are therefore reduced to 

T^G^\z,z'\zl,z2y^ = /o(^/?g")E^i2(-^i -^2)-^i?2i(^2-;^i)i'?Gf^(^,/ki*,^2)' 

k,l 

T2Gf\z,z'\zl,z2y^ = /o(^/lg")E ^12(^1 - Z2),'^i?2i(-^i - z2)i)G'^\z,z'\zlz2y 

k,l 

^ Note that there are typographical errors in the formula . 



(5.24) 



(5.25) 
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where zl — zi — ^w, and we use ( |2.28D and ( |3.1lD . Substituting ( ^.25D into ( ^.2lD we obtain 

P^\Z,Z'\Z,,Z2) = Cl^^A{u,)A{u2)Bra{u+)B _„j(u_). (5.26) 

Here Cm is a constant, and A{u) and Bm{u) are solutions to the foUowing difference equations: 

(5.27) 



Aiu) 



Brn{u) 



By solving these difference equations we obtain 
where 



2, (g2u4/n.i2^^4)^(^4^-4/n.^2^^4)^ 



(i2w4/«; g4)^(i2q2„-4/„. ^2^ ^4)^ ' 

i^{u) 5o(ug-"/2)5(^.-ig"/2), 5o(^i) 



(5.28) 



(q2+3«u-2. ^2^ ^4n)^(^2^-2+3ny-2. ^2^ ^2„^ ^4„)^ 



and 



where 



(g3«y2. ^2^ ^2«^ ^4n)_^(-^2^3«^2. ^2^ q2n ^ ^iny 
(^2^my2/„. ^2)^(^2^-m^-2/n. ^2)^ 



(g2tj"iu2/n. g2)_^(-^2^-my-2/«. ^2)^ ' 



(5.29) 



(g3"u-2; ^2^ ^2«^ ^2«-)^(-^2g3«y-2. ^2^ ^2«^ ^2n)^ 
/'^2+nT,2. +2 ^2n ^2n\ /'-f2^ — 2+n^,2. +2 ^2n ^2n\ 



Note that Bm{u) is essentially the same as G^"^\u) in which corresponds to the quantity (5.21) in 
the bulk theory. 



From ( 5.26|) we have 



Pi'\z,z' = 0\zi,Z2) 



Bi(u+)B_i(u_) 



pI,'\z,z' ^0\zi,Z2) C^'^ Bo{u+)Bo{u-) 



(5.30) 



Taking the low temperature limit g — > 0, we find that the ratio C^'^^Cq should be equal to lo^. We 
therefore obtain the boundary spontaneous polarization from (5.3C) and (5.29) by putting u+ = u_ = 1 

(5.31) 



(t2;t2)4^ (^2^.^2)2^ (52^-1.^4)2^- 

When n — 2 this expression coincides with the previous result obtained in We also emphasize that 
the boundary spontaneous polarization for the boundary Belavin model is exactly the square of that for 
the bulk Belavin model obtained in p2], up to a phase factor. 



6 Summary and discussion 

In this paper we have obtained two non-diagonal solutions of the reflection equation associated with 
Belavin's Z„-symmetric elliptic model. Unfortunately, our elliptic X-matrix is not connected with the 
diagonal boundary Boltzmann weig hts for the ^l^^l^-face model Ipl but the non-diagonal ones. It is 
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thus an open problem to obtain the X-niatrix correspondmg to the boundary Bohzniann weights given 
inH. 

On the basis of the boundary CTM bootstrap we have derived a set of difference equations for 
correlation functions of the boundary Belavin model. By solving the simplest difference equations, we 
have obtained the boundary spontaneous polarization of the boundary Belavin model. Our result is 
consistent with the one given in when n = 2. The boundary spontaneous polarization is equal to 
the square of the bulk spontaneous polarization ||2^ up to a phase factor. The same phenomena were 
observed in (|, |9| . 

In this paper we have shown that correlation functions of the boundary model satisfy the i?-matrix 
symmetry and the reflection properties, which are the boundary analogue of Smirnov's first two axioms 
[^8| . It may be interesting to construct integral formulae for correlation functions such that the integrand 
possesses the determinant structure as in Smirnov's integral [|8j. 

In [ p^ integral formulae for correlation functions of the boundary XY Z model by using bosonization 
of vertex operators |^ . In order to obtain the higher n generalization of , the construction of free 
field realization of the boundary Belavin model is required. It is a very hard but important work. 
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